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We study the radial Schrodinger equation for a particle of mass m in the field of a singular attrac- 
tive g 2 /r 4 potential with particular emphasis on the bound states problem. Using the regularization 
method of Beane et al, we solve analytically the corresponding "renormalization group flow" equa- 
tion. We find in agreement with previous studies that its solution exhibits a limit cycle behavior 
and has infinitely many branches. We show that a continuous choice for the solution corresponds 
to a given fixed number of bound states and to low energy phase shifts that vary continuously with 
energy. We study in detail the connection between this regularization method and a conventional 
method modifying the short range part of the potential with an infinitely repulsive hard core. We 
show that both methods yield bound states results in close agreement even though the regulariza- 
tion method of Beane et al. does not include explicitly any new scale in the problem. We further 
illustrate the use of the regularization method in the computation of electron bound states in the 
field of neutral polarizable molecules without dipole moment. We find the binding energy of s-wave 
polarization bound electrons in the field of Ceo molecules to be 17 meV for a scattering length 
corresponding to a hard core radius of the size of the molecule radius (~ 3.37 A). This result can be 
further compared with recent two-parameter fits using the Lennard- Jones potential yielding binding 
energies ranging from 3 to 25 meV. 

I. INTRODUCTION 

The renormalization of attractive singular potentials of the form 1/r™ with n > 2 was recently studied by Beane 
et al. We shall in the following refer to this renormalization method as the R- method. The purpose of this work 
is to analyze in more detail the case n — 4 with particular emphasis on the bound states problem. On the physical 
side, this potential describes the long range part of the polarizability potential in atomic and molecular systems and 
is relevant to the description of the long range proton-deuteron electromagnetic interaction. From a more formal 
viewpoint, it is of interest to study how the regularization method of Ref. |l| for the 1/r 4 potential compares with 
the results obtained for the 1 /r 2 potential 0, y] , and whether it agrees with previous renormalization schemes for 
the same interaction 0. 0] In this work, we follow Ref. in order to find an analytic form of the solution to the 
renormalization group flow equation. We then compute both the bound states spectrum and the low energy phase 
shifts arising from the renormalized potential and compare the R-method to a conventional method using a hard core 
radius for regularizing purposes. We then use the R-method to discuss the binding energy of s-wave polarization 
bound electrons in the field of neutral molecules with zero dipole moment. Our main results are the following. 

1. There are infinitely many solutions (3 n , [n — 1,2,3, . . .) to the renormalization group flow equation. Each (3 n 
exhibits a limit cycle behavior with, however, a period that depends on the cut-off radius R. Furthermore, (3 n 
takes the value nir in the limit R = 0. 

2. A continuous choice of solution, obtained by jumping from one branch to the next closest branch below with 
decreasing values of the cut-off radius, corresponds to a renormalization with a given fixed number of bound 
states. However only the energy level with the weakest binding energy is insensitive to the value of the cut-off 
radius. 

3. A numerical computation shows good agreements between the physical (in the sense discussed in Sec. Illlll bound 
states spectrum obtained with the R-method and the corresponding spectrum obtained in a conventional method 
parameterizing the scattering length with a hard core radius. When applied to the problem of bound electrons 
in the field of polarizable molecules without dipole moment, we find in particular that the R-method yields a 
binding energy of 17 meV for s-wave polarization bound electrons in the field of Cgo molecules. 
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Our paper is organized as follows. In Sec. [HI we present the R-method proposed in Ref. pj and we obtain the 
renormalization group flow equation for the singular attractive 1/r 4 potential. In Sec. lIIll we discuss the bound states 
spectrum of the regularized potential and discuss in Sec. IIVI its connection with a conventional method modifying 
the short range part of the potential by means of a hard core radius. In Sec. we discuss the application of the 
R-method to the calculation of electron binding energies in the field of polarizable neutral molecules with zero dipole 
moment. In Sec. lVIl we show that the low energy phase shifts are, as expected, insensitive to the cut-off radius. Some 
concluding remarks are reported in Sec. IVIII 



II. RENORMALIZATION METHOD 



In this paper, we follow the R-method proposed by Beane et al. pj to obtain analytically the renormalization group 
behavior of the coupling constant of the short-range attractive square-well used to regularize the singular attractive 
1/r 4 potential. 

We start with the s-wave reduced radial Schrodinger equation for one particle bound by a central potential V(r) 
(h = 2m=l): 



^2-V{r)-K 2 )r { r) 0. (I) 



with k = \J—E and where p|: 

V(r) = -^9(R-r)-^6(r-R) (a s ,a>0), (2) 

that is, the attractive (ai?) 2 /r 4 is cut-off at a short distance radius R by an attractive square well. As in Ref. pj, we 
first solve Eq. for the zero energy solution (k = 0) in order to find the corresponding wavefunction iPq(t). This 
solution is given by: 

ip (r) = Asm(a s ^) r<R (3) 

tpo(r) = B rcos [a— + <p\ r > R (4) 



where cj> is the zero energy phase p| and is given by 

tan^i = L/g (5) 

where L is the scattering length and g = aR. 

The usual matching condition of the wave function and its derivative at r = R then yields the renormalization 
group flow equation: 

a s cot a s = 1 + a tan(o + 4>). (6) 
We can solve analytically Eq. 10 as in Ref. to obtain: 

A) = ± iuJ ~ 1] ' cxp (- I' arg A (t)^ , « > 1 (7) 



.Wo t. 

(3 n = ±?i7rexp ( — f argO n (t) — ], — oo < uj < +oo, n = l,2, ... (8) 
Jo tJ 

where we denoted by (3 n the infinite set of solutions a s and we have: 

— = 1 + atanfa + (/>), (9) 

LU 

A (t) = X(t) + \wti-K, (10) 
X[t) = l + i^tln^f, (11) 

Q n (t) = A {tf +n 2 TT 2 uo 2 t 2 . (12) 
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FIG. 1: The running coupling constant f3 n as a function of a -1 = R/g for n = 1, 2, 3 and <j> = 1. The regions labelled £ = 1, 2, 3, 4 
are discussed in Sec. III. Quantities on both axes are dimensionless. 



The integer n is fixed on a given branch. The functions (3 n are given in Fig. ^for n = 1, 2,3 as a function of R. For 
computational ease, we chose <f> = 1. We only keep n > solutions as a; in formula (0 is unrestricted. It appears 
that the coupling constant of the square-well potential is a discontinuous function of R for a given n. A similar 
behavior was observed in the case of the renormalization of the singular 1 jr 2 potential. Note however two important 
differences: on the one hand, the period of oscillations, which changes with R, is no longer log-periodic 0,0,01; 
on the other hand, (3 n takes the value tvk in the limit R — > for all n. Indeed, the coupling constant, g = aR, of 
the 1/r 4 potential is fixed for a given physical system (for example g could be taken to be essentially the electrical 
polarizability, see Sec.[3>. Consequently, a varies with R and is infinite in the limit R — > 0. This leads to a vanishing 
u), see ©, and Q n (t) is then real yielding (3 n = nir (the sign of /?„, and thus of a s does not play a role, see JSJl). 

These results are consistent with the results of the study of the renormalization of long range attractive potentials 
in Ref. 4]. As already discussed in Ref. pj, one can also choose a s to be a continuous function of R. This implies 
jumping from one branch of the solution to the next one just below at the point of discontinuity as illustrated in 
Fig. n The respective merits of these two solutions were recently discussed in the literature in the case of the singular 
1/r 2 interaction HQ. 



III. BOUND STATES 



The renormalization method described in Sec. |H] only makes sense if the low energy observables (bound states and 
phase shifts for examples) are insensitive to the value of the cut-off radius R (for R small enough). The spectrum of 
the inverse square potential has been studied in detail Q, 0, an d it was shown that bound states with an energy 
above —1/R 2 are indeed insensitive to the value of the cut-off. 

As can be seen from Fig. ^ discontinuities of (3 n appear for singularities of u> (see (JSJl). For R large enough, the 
function |atan(a + 4>) \ is smaller than 1 and no singularity can appear for lu. The function (3 n is then continuous. This 
region is noted 1 in Fig.^ With the formula giving the number of bound states of the regularized potential obtained 
in Appendix [S] it is clear that in this region we have n bound states if the coupling constant of the square- well is 
obtained with (3 n (g — <f> —1). In the region noted 2 in Fig. ^ formula l|A5|l shows that the number of bound states 
has increased by one upon crossing a discontinuity of [3 n . Thus we must jump to a value (3 n ~i of a s in order to keep 
the number of bound states fixed. The various regions are obviously separated by the various discontinuities of (3 n . 
In general, in region £, the potential has n + i — 1 bound states for a s — (3 n (g = 4> = 1). Obviously, the number of 
bound states in the region i depends also on the value of g and <f> but this number always increases by one unit in the 
region i + 1 . 

With the analysis made above, it should be clear that a continuous choice of the solution a s (hence mixing several 
branches with different values of n) corresponds to a renormalization with a fixed number of bound states. This is 
illustrated in Fig. [3 where the evolution of the energy levels as a function of R is given when three bound states are 
present in the regularized potential. 

However, it is also clear that fixing the number of bound states in the potential yields a minimal value, R mm , for 
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FIG. 2: The bound states spectrum as a function of a 1 = R/g for a constant number of bound states fixed to three. Quantities 
on both axes are dimensionless. 



the cut-off radius R. For example, if we fix this number to 1 {g = <ft — 1), we must stay in the region 1 of Fig. ^ with 
a s computed with /3\. This leads to i? m in ~ 0.63 5. A crucial outcome of our calculations as shown in Fig. is the 
observation that the state with the weakest binding energy is insensitive to the value of R. The binding energy, as well 
as the mean square radius, of this state is also insensitive to the number of bound states present in the regularized 
potential. Since the states for which the binding energy varies with R have no physical meaning and since the binding 
energy and the mean square radius of the state with the weakest binding energy are insensitive to R and the number 
of bound states, there is no clear reason to choose a renormalization with a fixed number of bound states instead of 
a renormalization with a fixed branch since the latter does not introduce a minimal value for the cut-off radius R. 
Similar conclusions were drawn in the study of the singular attractive 1/r 2 potential 0, ■ We therefore conclude 
that the R-method used here yields bound states solutions independent of R and independent of the particular branch 
of the solution to the renormalization group equation. We shall thus refer to these renormalized solutions as physical, 
and dismiss the deeply bound, i?-dependent solutions as unphysical. 

Now, according to our numerical analysis, the binding energies Eb of these physical solutions for (f> — I are given 
by Eb = k 2 with k given by : 

3 k~0.83. (13) 

The form of formula l|13|) is counter-intuitive, as it implies that the renormalized binding energy increases with 
decreasing g. This behavior is nevertheless simple to understand since a decrease of g does not lead to an overall 
decrease of the regularized potential. The 1/r 4 part of the regularized potential, see is indeed less attractive but 
the square- well part of the potential, which depends also on g, see (JHJ), can be more attractive. The relation <|I3|) is 
consistent with a WKB analysis of the Schrodinger equation as shown in Appendix iBl 



IV. CONNECTION WITH HARD CORE POTENTIALS 



It is of interest to compare results from the R-method to those obtained with conventional methods where the 
polarizability potential — a p e 2 /(2r 4 ) is modified at the origin by means of some short range repulsion 0, Q , a p 
being the electric polarizability of the system. An especially simple form of such a modification is the hard core 
regularization which implies that the bound particle wavefunction must vanish at some hard core radius 1Z. The 
attractive feature of this model is that the corresponding scattering length L can be computed exactly to be : 




where ao is the Bohr radius. For a given value of L, Eq. (|14|) gives the corresponding value of 1Z for a given value of 
the electrical polarizability a p . On the other hand, the scattering length L determines the value of <f> from formula 
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TABLE I: Comparison between the binding energies Eb{= k 2 ) and the root mean square radii obtained with the R-method 
(with index R) and with the hard core regularization procedure (with index TV). The values of <f) an d TZ/g on the same line 
correspond to the same value of L from formulas JSJ and 1141 . 
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(J5J in the R-method. Comparison between Eqs. JSJl and l|14fl . and remembering that cj> is defined modulo ir, indicates 
that we have the correspondence 

K= (., + i/V-^ (15) 

with (f) e [0, tt] and s = 1, 2, 3, . . .. In Table [fl we show a comparison between the binding energies and root mean 
square radii obtained from both methods for the same values of the scattering length L, or equivalently for TZ and 4> 
related by (|15|) with s = 1. 

We see from TableQJthat the R-method and the hard core potential yield bound states values in excellent agreement. 
It is important to note, however, that this agreement holds for physical bound states, as defined above. This agreement 
also holds for any value of s. This integer is actually equal to the number of nodes of the wave function of the state 
with the weakest binding energy obtained with the hard core potential. However, there are significant differences 
between the two methods. The bound states wavefunction in the R-method has nodes corresponding to the number 
of unphysical (deeply bound) states and this number is actually arbitrary. The number of nodes of the hard core 
wavefunction is contrained by the value of 1Z which should correspond to some characteristic length of the system. 
The correspondence between the two methods leads to a further understanding of relation (TH^) . Indeed, TZ decreases 
with g in l|15fl leading to an increase of the binding energy. Moreover, since the number of nodes of the wave function 
obtained with the R-method varies with R (and the number of unphysical states), whereas the binding energy and 
mean square radius stay insensitive to this parameter, it is not easy to determine if the stable (physical) state is a 
ground state or an excited state. If the system considered has a given characteristic length, relations (|14f) and i|15|) 
can be used to determine s (since <fi is then known from (JjjJ) and then the position of the state in the spectrum. This 
is illustrated in the next section. 



V. POLARIZATION BOUND STATES 



In this section we apply the renormalization R-method to the computation of weakly bound electron states in the 
field of polarizable neutral molecules with zero dipole moment. The Cgo molecule is one of the few possible candidates 
with such properties 6] . In the conventional approach, bound states are computed by solving the Schrodinger equation 
with a 2-parameter Lennard- Jones potential {7jf. Taking the electrical polarizability a p of Cqo to be 558 Oq leads to 
an electron binding energy between 3 and 25 meV according to the value of the parameters describing the short 
range part of the interaction. We find a physical bound state at 17 meV with <fi = 1 corresponding to a scattering 
length L = gj tan^ with g = V ma p e 2 = aoV558. In a conventional hard core analysis, a p e 2 /2 = 279 e 2 ap, the same 
binding energy is obtained with an effective radius of 6.37 ao (~ 3.37 A) for the Ceo molecule. This radius is obtained 
with the relation (|15|l with 8=1. Larger values of s yield the same binding energy but the effective radius has no 
longer a physical interpretation. Moreover, this weakly bound state is certainly a ground state \J\. The value of the 
effective radius found with l)15ll is close to the experimental value of the mean radius of the Ceo molecule 3.55 A [To|. 
If this experimental value is used as effective radius, we find with (|15|l <fi — 1.192 and the binding energy is then equal 
to 6 meV (the relation ljl"3*|l becomes for this value of 0, gn ~ 0.50). Note that, as already clear from the results of 
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Table H] the binding energy depends strongly on the value of </>. Thus the scattering length, which determines the 
value of cj> (see iJSJ) should be known with enough precision to allow definite predictions and comparisons with other 
regularization methods. 



VI. LOW ENERGY PHASE SHIFTS 



Another test of the method described in Sec. [n] is the computation of the low energy phase shift. In Fig. [3] we 
plot the s-wave phase shift So as a function of gk (k = y/E > 0) for different values of R. One can see again that 
So is insensitive to the value of R for gk < 1, even though the number of bound states (physical and nonphysical, as 
defined above) varies from 1 to 4. From the viewpoint of the renormalization method, only one (weakly bound) state 
is physical and leads to Sq — tt for a vanishing energy. 



VII. CONCLUSIONS 



In this paper we studied, with due emphasis on the bound states problem, the renormalization of the singular 1/r 4 
potential using the method of Beane et al lj. We found, in agreement with previous works, that the solution (3 n to 
renormalization group flow equation exhibits a limit cycle behavior and has infinitely many branches. We discussed 
the form of the bound states spectrum as a function of f3 n and discussed the connection between the R-method and a 
conventional method using a hard core radius to modify the short range part of the interaction. We then applied the 
R-method to calculate the energy of polarization bound electrons to a neutral polarizable molecule without dipole 
moment. When applied to the Cgo molecule, we found that only very accurate values for the scattering length could 
discriminate between results from the R-method and those obtained from a polarization potential with hard core. 



APPENDIX A: NUMBER OF BOUND STATES OF THE RENORMALIZED POTENTIAL 



It is well known that the number of bound states in a central potential is equal to the number of zeroes of the zero 
energy wave function in the interval < r < oo (see for example |TljV Equivalently, we count the number of extrema 
of this wave function. Thus we have to count the number of zeroes of the derivative of the zero energy wave function 

V>o(r) ~ cos(a s -^) for r < R, (Al) 

V»o(r) ~ cos +<?!>) +^ sin (2 + ^) for r > R (A2) 

The number of extrema, N± , of the zero energy wave function in the interval < r < R is simply given 

*-{{? + !}}. (A3) 
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where {{2}} is the integer part of x. To find the number of extrema, N2, of the zero energy wave function in the 
interval R < r < 00, we count the number of zeroes of the expression cos(.t + <j)) + x s\xi(x + </>) for x < a. Equivalcntly 
we search for the number of solutions of the equation 



xt&n(x 



-1 x < a. 



(A4) 



The solutions are the intersection of a tangent function with a hyperbola. The number of solutions of the equation 
(|A4(I is given by 



N 2 = U- a 



arctan — 
a 



The total number is thus given by N — Ni + N2 with N% and N2 defined by (|A3jl and (|A5|I respectively. 

APPENDIX B: WKB ANALYSIS 

The general formula giving the energy spectrum of a central potential in the WKB approximation is: 



(A5) 



drsjE - V(r) = (n - 1/2)tt, 



(Bl) 



where n = 1, 2, 3, ... and r± are the solution to E = V(r±). We then get from Eq. @: 



R 

Performing the integrations in Eq. (|B2|I and writing x — R/r+, we find: 



v/£i? 2 + (a s ) 2 + ^ 



2r(5/4) 2 



= (n - 1/2)tt, 



(B2) 



(B3) 



where a, b) is the Beta function. Assuming that E remains finite as R — > and using the formulas for x — > 0: 



1 x s 



4x d 



V-l + l/^ 4 ^--y and i?(x 4 , 3/4,1/2)- — , 



we eventually get, keeping the leading term as R — ► 0: 



r(5/4) 



r(3/4) 



a s + ^ - 7r(n - 1/2) 
it 



(B4) 



(B5) 



From Eq. (| A3|> and Eq. (|A5|) . we find that n7r ~ a s + g/R + <j> in the limit R — > 0, so that we finally get from Eq. (|B5 



7K = 5TV--B 



r(5/4) 



r(3/4) 



1/2] 2 



(B6) 



Equation l|B6|l has the same functional form as Eq. (|13|l . It does not of course yield the same numerical value for the 
binding energy. 
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